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( U)  This  report  addresses  the  problem  of  correlating  explosive  reaction  in  a  warhead  penetrating  a  harrier 
to  the  stress  environment  imposed  on  it.  A  missile  fired  against  a  ship  and  a  missile  fired  against  a  concrete 
structure  represent  the  penetrating  impacts  of  special  interest  The  emphasis  of  this  work  is  on  the  response  of  the 
explosive.  A  standard  warhead  design  was  selected  because  a  predictive  model  was  developed  from  theories  of 
shock  generation,  transmission  and  reflection  in  metals,  and  existing  concepts  of  initiation  of  reactions  varying 
from  burning  to  detonation.  The  model  uses  equation  of  state  data  for  the  warhead  case,  barrier  and  explosive,  and 
wedge  test  data  and  underwater  sensitivity  test  data  to  predict  critical  combinations  of  impact  speed  and  barrier 
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INTRODUCTION 

The  defeat  of  targets  that  have  a  protective  barrier  enclosing  vulnerable  components 
requires  that,  in  order  to  achieve  optimum  effectiveness,  a  missile  warhead  fully  penetrate 
the  barrier  prior  to  detonation.  The  warhead  must  not  only  pass  through  the  barrier  but  also 
must  withstand  the  dynamic  loading  imposed  on  it  without  losing  the  functional  integrity  of 
essential  components.  In  particular,  the  explosive  charge  must  survive  the  dynamic  loading 
without  significant  reaction.  Two  examples  are  of  particular  interest:  (1)A  missile  fired 
against  a  ship  target  achieves  the  greatest  damage  to  the  ship  if  the  missile  perforates  the 
hull  and  penetrates  into  the  ship’s  interior  without  premature  initiation  of  the  explosive 
charge  so  as  to  take  advantage  of  the  enhancement  of  damage  from  confinement.  (2)  A  missile 
fired  against  a  concrete  structure  needs  to  penetrate  into  the  barrier  material  before 
detonation  in  order  to  achieve  confinement  and  to  avoid  attenuation  of  damage  by  the  barrier 
itself. 

This  report  addresses  the  problem  of  correlating  explosive  reaction  to  the  stress 
environment  imposed  on  it  by  a  penetrating  impact.  The  process  of  correlation  has  two  parts 
The  first  part  looks  at  interaction  between  the  structural  components  of  the  warhead  and  the 
barrier.  Specifically,  the  generation  of  stresses  bv  impact  and  the  subsequent  dynamic 
loading  of  the  explosive.  The  second  part  concerns  the  response  of  the  explosive  to  dynamic 
loading  The  emphasis  of  this  report  will  be  on  the  second  part.  Therefore,  the  investigation 
has  minimized  the  variations  in  warhead  configuration  in  order  to  achieve  maximum  clarity 
in  the  characterization  of  explosive  response 

To  say  that  a  warhead  design  achieves  clarity  in  the  characterization  of  explosive 
response  to  dynamic  loading  means  that  in  experiments  (actual  or  hypothetical)  the  design 
develops  all  of  the  relevant  and  representative  aspects  necessary  for  the  desired  correlation, 
and  the  specific  details  of  its  configuration  and  material  properties  are  directly  related  to  the 
dynamic  loading  The  section  entitled  Background  and  Approach  discusses  the  dynamics  of 
impact  and  explosive  sensitivity  and  shows  that  structurally  robust  flat-ended  warheads, 
striking  steel  plates  and  concrete  slabs  at  normal  incidence,  contribute  to  this  kind  of  clarity. 

The  section  entitled  Predictive  Model  describes  models  for  predicting  premature 
reactions  in  the  chosen  warhead  designs  for  penetrating  impacts  against  steel  plates  and 
concrete  slabs  These  models  integrate  the  theories  and  data  discussed  in  Background  and 
Approach  into  a  series  of  comparatively  simple  analytical  procedures  by  which  one  can 
predict  the  conditions  for  premature  reactions  given  the  impact  conditions  and  warhead  and 
target  parameters. 
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BACKGROUND  AND  APPROACH 


DYNAMIC  LOADING  OF  PENETRATOR  AND 
ITS  EXPLOSIVE  COMPONENTS 

The  dynamic  loading  of  the  explosive  components  of  a  penetrating  warhead  involves  a 
sequence  of  events  that  (1)  begin  with  the  generation  of  stresses  at  the  interface  between 
penetrator  and  barrier  (2)  followed  by  the  propagation  of  a  distribution  of  stresses  throughout 
the  barrier  and  the  penetrator,  and  (3)  in  particular,  the  application  of  dynamic  loading  to  the 
explosive-filled  region.  Mathematical  formulations  for  the  generation  of  stresses  at  the 
interface  between  projectile  and  target  were  part  of  the  earliest  work  in  terminal  ballistics 
(References  1  through  4).  In  addition,  more  recent  formulations  have  had  at  least  limited 
success  in  predicting  the  motion  of  rigid  projectiles  from  specified  behavior  of  deforming 
media  (References  5  through  7).  With  advances  in  computer  technology,  finite  element  and 
finite  difference  formulations  of  impact  problems  in  terms  of  continuum  mechanical 
principles  have  taken  advantage  of  high-speed  computers  to  do  the  horrendous  bookkeeping 
involved  in  very  fundamental  analyses  of  impact  (References  8  through  10).  These  analyses 
pose  the  impact  problems  as  boundary  value  problems  for  the  partial  differential  equations 
that  express  the  basic  conservation  laws  of  physics  and  a  constitutive  equation  that  defines 
material  properties.  This  type  of  analysis  gives  solutions  for  virtually  arbitrary 
configurations  and  material  properties,  and  the  solutions  include  the  details  of  the 
development  of  internal  stresses  and  strains.  These  predictions  of  impact  phenomena  agree 
with  the  trends  of  experimental  observations  and  provide  a  useful  basis  for  planning  and 
interpreting  experiments.  Solutions  by  this  method  affirm  the  early  dominance  of  wave 
propagation  in  the  distribution  of  stresses  and  the  simplification  of  phenomena  for  minimal 
deformation  of  the  penetrator.  Such  solutions  also  substantiate  the  considerable  advantage 
derived  from  the  use  of  a  flat-ended  penetrator. 

A  flat-ended  penetrator  striking  a  plate  at  normal  incidence  simultaneously  produces 
a  single  mechanical  state  over  the  entire  surface  of  contact.  All  other  shapes,  such  as  conical 
or  ogival,  develop  a  contact  surface  at  a  finite  rate  of  the  order  of  magnitude  of  the  penetration 
rate.  The  interactions  between  the  surfaces  have  varying  degrees  of  development  along  the 
surface  so  that,  for  example,  the  stresses  from  the  tip  of  the  pointed  penetrator  have 
propagated  well  into  the  warhead  while  the  stresses  from  the  most  recent  contact  lie  close  to 
the  penetrator. 

On  the  other  hand,  the  planar  symmetry  of  the  contact  surface  between  the  flat-ended 
penetrator  and  plate  ensures  initial  states  of  stress  and  strain  that  depend  on  only  one  space 
dimension,  and  the  propagation  of  this  initial  state  becomes  an  exercise  in  one-dimensional 
wave  propagation.  The  major  analytical  task  becomes  that  of  representing  the  material 
responses  to  the  collision  of  surface.  Unfortunately,  this  initial  state  of  great  simplicity  soon 
fades  with  the  propagation  of  relief  effects  from  the  lateral  free  surfaces.  A  very  crude  and 
simple  calculation  of  the  rate  of  traversal  of  the  barrier  and  nose  plate  shows  that  if  the  ratio 
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where 

hb  =  the  thickness  of  the  barrier 
hn  =  the  thickness  of  the  warhead  nose 
D  =  the  diameter  of  the  projectile 

is  applied,  then  parts  of  the  explosive  experience  only  this  simple  form  of  dynamic  loading  in 
the  early  part  of  the  impact. 
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Clearly,  any  system,  no  matter  how  esoteric,  that  exhibits  such  simplifications  of 
behavior  has  considerable  merit  as  the  system  of  choice  in  an  investigation  of  the  sensitivity 
of  internal  explosive  components  to  dynamic  loading  from  impact.  It  turns  out  that  recent 
designs  of  penetrating  ordnance  have  chosen  the  flat  ended  configuration  for  very  practical 
reasons  unrelated  to  explosive  sensitivity,  e  g  ,  to  achieve  yaw  stability  in  the  penetration  of 
extended  targets  (Reference  11).  Experimental  studies  of  the  stability  of  penetrators  have 
shown  that  instability  becomes  a  minimum  for  the  flat  ended  nose  shape  and  other 
investigations  show  that  pointed  shapes  worsen  instability  bv  encouraging  cavitational 
processes  that  put  the  center  of  pressure  forward  and  thus  develop  an  unstable  condition  The 
corner  at  the  edge  of  the  flat-ended  penetrator  has  a  strong  tendency  to  increase  the  efficiency 
in  the  penetration  of  thin  plates  and  it  also  greatly  reduces  the  tendency  to  ricochet 
(References  12  and  13).  Thus,  the  fiat-ended  nose  shape,  which  promises  such  significant 
simplification  of  the  analysis  of  internal  loading,  also  has  direct  relevance  to  existing 
ordnance.  For  these  reasons,  the  model  described  here  refers  exclusively  to  a  fiat-ended 
cylindrical  penetrator  with  a  flat-ended  cylindrical  explosive-filled  cavity.  The  model’s 
predictions  come  from  the  application  of  plane  shock  wave  theory  to  the  initial  stages  of  the 
dynamic  loading  developed  in  the  penetrating  impact. 


PLANE  SHOCK  WAVES 

For  many  years,  experiments  using  plane  shock  propagation  have  provided  most  of 
the  information  on  the  behavior  of  materials  at  very  high  pressure  and  very  high  strain  rate. 
Explosive  plane  wave  generators  (Reference  14)  or  the  impact  of  cylinders  against  plate 
(Reference  15)  have  developed  plane  shock  waves  in  materials  of  interest  Measurements  of 
the  propagation  rate  and  particle  speeds  of  these  shocks  provide  the  data  needed  to 
characterize  the  materials.  The  conventional  theory  of  shocks  in  fluids  provides  the  basis  for 
the  interpretation  of  these  experimental  measurements  The  equations  for  fluids  apply  to  the 
solid  target  and  warhead  case  on  the  premise  that  the  stresses  in  the  material  exceed  its  yield 
strength  so  that  the  material  has  lost  resistance  to  shear  but  not  to  compression.  The 
behavior  of  the  material  then  strongly  resembles  that  of  a  fluid 
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The  governing  equations  for  one-dimensional  shocks  come  from  the  application  of  the 
principles  of  the  conservation  of  mass,  momentum,  and  energy  to  the  flow  through  the  shock 
front.  These  conservation  equations  take  the  form  given  below  (Reference  16). 
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E-E  =  -(p+P  )(  — 
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where 

U  =  shock  speed 
u  =  particle  speed 
P  =  pressure  behind  a  shock  front 
P0  =  pressure  ahead  of  a  shock  front 
E  =  internal  energy  behind  shock  front 
E0  =  internal  energy  ahead  of  shock  front 
p  =  density  behind  shock  front 
p„  =  density  ahead  of  shock  front 


(3) 


Equations  1  and  2  suffice  to  establish  the  dependence  of  the  shock  speed  and  the  particle  speed 
(u)  of  the  parameters  U  and  p.  Equation  3,  called  the  Hugoniot  equation,  expresses  the 
unique  properties  of  specific  material.  The  data  needed  to  determine  a  particular  Hugoniot 
equation  may  come  from  either  the  equation  of  state  of  the  material  or  from  measurements  of 
U  and  u  in  plane  shock  experiments.  For  many  materials,  the  shock  and  particle  speeds  have 
a  linear  relationship 

U  =  a  +  bu  (4) 


The  combination  of  Equations  2  and  4  determine  the  dependence  of  P  on  u.  The  equation  uses 
the  same  information  as  the  Hugoniot  equation  and  comprises  an  equivalent  representation 
of  the  behavior  of  a  particular  material.  The  equation  in  the  form 


P  =  p(a  4-  bu)u  (5) 

has  the  common  name  "the  direct  Hugoniot”  to  distinguish  it  from  "the  reflection  Hugoniot,” 
which  will  be  discussed  in  a  later  section. 


Shock  Generation 

Two  bodies  that  collide  on  plane  surfaces  with  an  impact  speed,  V,  abruptly  develop  a 
common  pressure,  P,  and  a  common  surface  motion,  u,  over  the  surface  of  contact.  The 
pressure  and  the  surface  speed  belong  to  wave  motions  that  propagate  into  each  of  the 
colliding  bodies.  These  wave  motions  serve  to  accommodate  the  two  bodies  to  the  differences 
in  motion  after  contact  occurs.  Due  to  the  abruptness  and  intensity  of  the  processes,  the 
waves  have  the  form  of  a  shock.  The  shocks  have  a  common  pressure,  but  each  has  a  particle 
speed  that  depends  on  the  Hugoniot,  Equation  3,  for  the  particular  material 
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The  particle  speeds  also  account  for  the  closure  speed  so  that 


V  =  Uj  + 


(9) 


therefore 


P  =  P1tf1u1  =  p2l/2(V-u1)  (10) 

Figure  1  graphically  presents  this  equation.  The  two  curves  represent  the  direct 
Hugoniot  of  material  one  and  the  reflection  Hugoniot  of  material  two.  The  latter  curve 
represents  the  conditions  given  by 


P  =  P2^2('/  ~  "i*  =  P2a2iV  ~  “l*  +  P2h2(V  ~  “l)2  (11> 

The  shock  motion  in  this  second  body  has  the  opposite  direction  from  that  of  the  first  body  so 
that  it  has  a  mirror  image  dependence  on  the  particle  speed,  U  i ,  with  the  intersection  of  the  U 
axis  at  Ui  =  V.  Such  curves  represent  the  same  information  as  the  Hugoniot  with  the  form  of 
Equation  3,  but  it  becomes  clear  that  from  equations  such  as  Equation  11,  the  representation 
depends  on  the  reference  system  used  for  the  observation  of  the  particle  speed.  This  kind  of 
situation,  in  which  it  becomes  convenient  to  compare  particle  speed  of  different  directions 
with  respect  to  a  single  reference  frame,  occurs  frequently  in  problems  of  shock  reflection. 


FIGURE  1.  Examples  of  a  Direct  Hugoniot  and  a 
Reflection  Hugoniot. 
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Clearly,  if  the  nose  and  barrier  materials  have  identical  Hugoniot  representations  (of 
the  form  of  Equation  3),  then  the  two  curves  of  Figure  1  have  a  mirror  image  relationship  and 
the  particle  speed  becomes 


and  the  pressure 
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p  =  pi,;i7 


(12) 


(13) 


Shock  Attenuation  and  Wave  Shape 

At  very  high  stresses,  metals  behave  according  to  the  laws  of  plasticity  rather  than 
the  laws  of  elasticity  Shear  components  of  stress  have  limitations  imposed  by  the  effects  of 
permanent  deformation  that  have  the  properties  of  flow  and  give  the  material  behavior  a 
resemblance  to  fluid  flow  The  resemblance  to  fluid  flow  increases  with  the  intensity  of  the 
applied  loads  and  justifies  the  use  of  fluid  theory  of  shocks  for  impacts  greater  than  a  few 
hundred  meters  per  second.  Equations  1  through  3  provide  the  basis  for  determining  the 
pressure  and  particle  speed  of  shock  fronts  generated  by  impact  according  to  Equations  5 
through  9.  These  also  form  the  basis  for  the  next  section's  discussion  of  the  phenomena  at 
interfaces  between  materials  of  different  properties  and  interactions  among  shocks.  These 
predictions  all  agree  with  carefully  made  observations;  however,  the  same  kind  of  careful 
observations  show  that  other  peculiarly  solid  phenomena  enter  the  process  and  result  in  a 
multiwave  system  involving  more  than  one  propagation  rate  and  the  resultant  changes  of 
wave  shape.  Figure  2  shows  an  experimental  measurement  of  the  pressure-time  profile  made 


FIGURE  2.  Pressure-Time  Profile  Measured  at  the  Rear 
Surface  of  an  Explosively-Loaded  Mild  Steel  Plate. 
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at  the  rear  surface  of  a  plate  of  mild  steel  explosively  loaded  on  the  front  surface 
(Reference  17).  This  wave  has  a  typical  multiwave  structure.  The  two  step-like  features 
propagate  at  different  rates.  The  wave  of  lower  pressure  propagates  at  the  speed  of  elastic 
disturbances  and  has  magnitude  approximately  equal  to  the  dynamic  yield  value  of  the 
material  and,  thus,  corresponds  to  the  residual  elastic  effects  in  the  plastically  deforming 
body  (Reference  18).  Phase  changes  in  the  solid  will  also  produce  additional  structure  in  the 
forward  part  of  the  wave  shape  (Reference  19).  In  summary,  shock  phenomena  in  metals 
follow  fluid  dynamic  principles,  but  the  full-wave  system  exhibits  the  effects  of  residual 
elasticity  and  phase  changes  characteristic  of  the  solid  material. 

The  development  of  the  multiwave  structure  just  described  and  the  dissipative 
processes  that  occur  in  a  shock  result  in  the  attenuation  of  the  wave  so  that  as  it  propagates, 
its  shape  and  its  amplitude  change  The  original  pressure,  P,  developed  by  impact  decays  to  a 
lower  value,  P.  An  expression  for  the  attenuation  at  a  distance  of  propagation,  S.  based  on 
simple  proportionality  of  the  rate  of  loss  of  intensity  to  the  instantaneous  intensity  has  the 
form 

P  ,=Pe-aS  <14) 

Experience  shows  that  this  simple  form  adequately  describes  the  net  attenuation  of  many 
materials.  In  this  simplified  description,  the  one  constant,  a,  covers  the  net  effect  of  both  the 
development  of  the  multiwave  structure  and  dissipation  due  to  the  irreversible  processes 
active  in  the  shock,  such  as  heat  conduction  and  internal  viscosity. 


Shock  Reflection  and  Transmission 

A  discontinuity  exists  in  material  properties  at  an  interface  between  two  different 
material0  A  shock  arriving  at  such  an  interface  encounters  a  discontinuity  in  the  Hugoniot 
that  determines  the  relation  between  pressure  and  particle  speed.  Despite  the  discontinuous 
change  in  material  properties,  the  dynamic  state  must  remain  continuous  except  at  a  shock 
front.  This  can  occur  by  generation  of  a  reflected  shock  in  the  material  in  which  the  incident 
shock  approached  the  interface  and  a  transmitted  shock  in  the  second  material.  These  shocks 
establish  a  new  common  pressure  at  the  interface  and  particle  speeds  in  each  body  consistent 
with  a  common  motion  at  the  interface. 

The  Hugoniot  for  the  state  of  the  first  material  after  generation  of  the  reflected  shock 
and  in  the  pressure-particle-speed  form  (and  in  the  frame  of  reference  for  the  incident  shock) 
becomes  the  mirror  image  of  the  direct  Hugoniot  that  also  passes  through  the  pressure  and 
particle  speed  of  the  incident  shock  (Figure  3).  This  now  represents  the  state  achieved  in  the 
first  material  under  the  combined  incident  and  reflected  waves  and  with  a  surface  motion  of 
the  interface,  U2,  (in  the  original  frame  of  reference  ).  The  equation  in  u>  for  a  common 
pressure  at  the  interface  becomes 


P  =  P2U2U2  =  P\U\i2ui  ~  u 


1  1 


Pt0  ,(2w  , 


It,,)  + 
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FIGURE  3.  Hugoniot  Curves  Representing  the  Shock  Generation, 
Reflection,  and  Transmission  for  an  Impact  at  Two  Plane  Surfaces. 


\s  here 


f/j  =  a(  +  b  i2u  —  «.,) 

U •?  =  a2  +  b2U2 

The  expression  on  the  far  right  side  of  Equation  15  represents  the  reflection  Hugoniot,  i.e.,it 
describes  the  state  of  the  first  material  that  satisfied  the  Hugoniot  of  the  first  material  (in  its 
P  -  p  form  of  Equation  3),  but  corresponds  to  surface  motions,  u-2,  under  these  new  conditions 
of  combined  incident  and  reflected  shock.  Under  these  new  conditions,  the  particle  speed 
contributed  from  the  reflected  shock  corresponds  to  particle  speed  with  the  opposite  effect  on 
the  pressure  because  reflection  has  reversed  the  direction  of  propagation  of  the  reflected 
shock  The  intersection  of  this  reflected  Hugoniot  and  the  direct  Hugoniot  of  the  second 
material  established  the  particular  motion  of  the  interface  compatible  with  the  discontinuity 
of  material  properties. 

This  kind  of  process  occurs  at  the  interface  between  the  nose  plate  material  and  the 
explosive,  and  the  above  principles  of  "impedance  matching”  determine  the  amplitude  of  the 
shock  that  enters  the  explosive  If  the  first  material  has  a  free  surface  at  which  no  significant 
pressure  can  develop,  the  Hugoniot  for  the  "second  material”  collapses  to  the  u-axis,  the 
particle  speed  becomes  twice  the  incident  particle  speed,  and  the  pressure  after  reflection 
becomes  zero.  A  reflection  of  this  kind  cancels  the  pressure  developed  at  impact  and  puts  a 
finite  duration  on  the  stress  wave  due  to  impact. 
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EXPLOSIVE  SENSITIVITY  TO  SHOCK  LOADING 

An  explosive  responds  to  shock  wave  loading  in  a  variety  of  ways.  (1 )  The  explosive 
may  detonate  almost  instantaneously.  (2)  A  low-level  chemical  reaction  may  occur  and 
develop  into  a  detonation.  (3)  A  low-level  reaction  may  occur  then  either  dies  out  or  builds  in 
intensity  to  a  destructive  level,  but  not  to  a  true  detonation.  Many  other  observations  may 
reflect  the  importance  of  some  particular  characteristic  of  the  warhead,  such  as  the  degree  of 
confinement  of  the  region  of  reaction.  The  particular  response  does  depend  on  one  or  more  of 
the  parameters  that  characterize  the  shock  wave:  P,  U,  u,  or  T.  At  high  shock  pressure 
greater  than  a  critical  pressure,  Pj,  the  shock  becomes  a  detonation  wave  in  times  too  short  to 
resolve  with  current  measurement  techniques. 

Above  a  significant  lower  threshhold,  Pc,  the  shock  and  the  chemical  reaction  that  it 
induces  develop  as  distinct  but  interrelated  phenomena.  The  shock  builds  up  due  to  the 
energy  fed  into  it  by  the  chemical  reaction,  and  the  chemical  reaction  builds  up  due  to  the 
greater  intensity  of  the  shock.  The  reaction  zone  chases  the  shock  as  each  phenomena 
augments  the  other.  The  process  of  mutual  augmentation  continues  until  it  reaches  the  well- 
known  limit  at  the  conditions  for  a  steady  detonation  wave  This  occurs  after  the  shock  and 
reaction  zone  have  traveled  a  distance  that  depends  on  both  the  initial  shock  conditions  and 
the  explosive,  the  distance  becoming  progressively  longer  for  lower  shock  intensities.  Beyond 
a  distance  of  the  order  of  20-30  mm,  both  the  shock  and  the  reaction  zone  die  out.  If  any 
sustained  reaction  occurs,  it  has  a  different  nature.  The  kind  of  buildup  of  reaction  just 
described  has  the  designation  shock-to-detonation  transfer  (SI)T)  ( References  20,  2 1 .  and  22) 

At  still  lower  pressures,  less  than  a  critical  pressure  (Pc),  the  shock  and  chemical 
reaction  become  completely  distinct  phenmena  once  the  shock  has  initiated  the  chemical 
reaction.  Experimental  studies  in  this  regime,  called  the  low-amplitude  long  duration  shock 
regime  (LALDS),  indicate  that  the  impulse  in  the  shock  wave  correlates  to  the  initiation  of 
sustained  burning.  Hence,  the  shock  pressure,  the  pulse  duration,  and  the  wave  form  all 
enter  into  the  critical  conditions  for  the  initiation  of  burning  (Reference  23). 

Shock-to- Detonation  Transfer 

Investigations  in  the  SDT  regime  have  developed  a  criterion  for  the  transition  to 
detonation  that  takes  the  form  of  a  critical  energyfluence  criterion  (Reference  24)  given  by 

p2y  _  ^  (16) 

where  the  factor  T  represents  the  induction  time  for  the  reaction,  i.e. ,  the  times  between  the 
arrival  of  the  pressure  and  temperatures  conditions  of  the  shock  front  at  a  given  point  in  the 
explosive  and  the  manifestation  of  a  significant  level  of  reaction  at  that  same  point  Some 
have  defined  the  level  of  reaction  as  the  release  of  1%  of  the  total  available  chemical  energy. 

The  minimum  pressure  for  which  this  kind  of  reaction  occurs  corresponds,  by 
Equation  16,  to  very  short  times.  This  virtually  guarantees  that  the  duration  of  the  pressure 
will  meet  the  time  requirements  for  any  meaningful  barrier  thickness.  1'hus,  although  the 
critical  criterion  for  initiation  does  depend  on  time,  the  brevity  of  the  times  for  all  practical 
purposes  makes  the  initiation  dependent  only  on  exceeding  a  critical  pressure,  Pr 
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Reaction  From  Large-Amplitude 
Long-Duration  Shocks 

Below  the  SDT  range  of  pressures  is  the  large-amplitude  long-duration  shock 
(LALDS)  regime.  In  the  LALDS  regime,  an  equation  similar  to  Equation  16  serves  to 
separate  conditions  of  successful  from  unsuccessful  ignition  (Reference  15): 


The  difference  in  form  has  no  ready  explanation  as  a  phenomenon  except  that  it  correspond  to 
broad  changes  in  the  physical  processes  responsible  for  ignition. 

As  mentioned  above,  reactions  in  the  LALDS  regime  show  no  sustained  coupling 
between  the  initiating  shock  and  the  reacting  region  of  the  explosive.  The  reaction  exhibits 
accelerated  growth  but  at  a  far  slower  rate  than  for  SDT,  and  the  final  state  of  reaction 
remains  separated  from  the  initiating  shock  and  need  never  have  the  properties  of  a 
detonation.  Whereas,  SDT  appears  to  result  from  conditions  throughout  the  region  on  or  near 
the  shock  front  that  favor  partial  reaction  of  the  explosive  material,  LALDS  involves  highly 
localized  responses  at  very  small  voids,  inclusions,  and  perhaps  other  inhomogeneities.  A 
considerable  amount  of  research  on  the  mechanisms  of  responses  to  shock  in  this  regime  has 
established  the  involvement  of  minute  structural  or  compositional  flaws  and  has  provided 
models  for  coupling  of  energy  from  the  shock  into  the  explosive  material  (Reference  26).  For 
the  purposes  of  this  study,  one  needs  a  quantification  of  the  submacroscopic  processes  in  terms 
of  macroscopic  material  properties.  For  example,  the  number  and  size  of  voids  correspond  to 
changes  in  density.  Thus,  the  sensitivity  of  an  explosive,  as  measured  by  the  critical  pressure 
(Pc),  does  in  fact  increase  dramatically  with  decreases  of  density  from  the  theoretical 
maximum  density.  Other  characteristics  of  the  internal  structure  involve  relatively  complex 
measurements  that  for  most  materials  do  not  exist.  In  general,  the  models  for  explosive 
response  to  shock  loading  prove  inadequate  to  make  predictions  in  detail,  and  would  require 
information  from  expensive  tests. 

Experiments  such  as  the  underwater  sensitivity  test  (Reference  23)  establish  the 
dependence  of  ignition  on  pressure  duration.  Apparently  the  lower  pressures  cannot  produce 
reactions  that  have  sufficient  intensity  to  propagate  forward  with  a  net  contribution  to  the 
shock  strength.  However,  given  an  apropriate  duration,  these  pressures  can  produce 
reactions  of  sufficient  intensity  to  result  in  growth  and  coalescence  of  the  system  of  localized 
reaction  into  a  self-sustaining  burn. 


CHOICE  OF  WARHEAD-BARRIER  SYSTEM 

The  susceptibility  of  explosives  to  premature  reactions  during  a  penetrating  impact 
depends  on  the  dynamic  loading  applied  by  the  warhead  to  the  explosive  and  on  the  particular 
response  of  the  explosive  to  this  loading  The  interaction  of  the  warhead  body  and  the  barrier 
determines  the  dynamic  loading  given  to  the  explosive.  This  study  focuses  on  (1)  how  an 
explosive  reacts  to  typical  loading  and  (2)  the  comparison  of  various  explosives.  Therefore, 
the  warhead  configuration  should  remain  as  constant  as  possible.  Furthermore,  the  chosen 
warhead  configurations  should  adequately  represent  real  configurations  and  allow  the  most 
simple  and  direct  correlation  of  dynamic  loading  to  the  warhead  and  to  impact  parameters. 
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The  preceding  sections  have  shown  the  feasibility  of  such  correlations  using  planar 
symmetry,  e.g.,  a  flat-ended  warhead  impacting  a  plate  at  normal  incidence,  which  represents 
the  conditions  encountered  by  in-service  warheads. 

The  considerations  of  plane  shock  phenomena  and  explosive  sensitivity  suggest  the 
following  limitations  on  the  impact  system  in  order  to  make  the  analyses  of  dynamic  loading 
meaningful  and  tractable: 

1.  The  warhead  has  a  flat  nose  configuration  that  presents  a  flat  surface  for  contact 
with  the  barrier  and  has  a  flat-ended  explosve-filled  chamber  with  the  warhead- 
explosive  interface  parallel  to  the  plane  of  the  nose. 

2.  The  warhead  has  material  strength  and  wall  and  nose  plate  dimensions  large 
enough  to  ensure  minimal  deformation  of  the  warhead  case. 

3.  The  target  barrier  consists  of  either  a  plate  of  steel  or  a  slab  of  concrete. 

4.  A  steel  barrier  has  a  Hugoniot  identical  to  that  of  the  warhead  nose  plate. 

ASSUMPTIONS  AND  SIMPLIFICATION 

Given  the  above  warhead,  the  model  uses  certain  assumptions  to  simplify  and 
expedite  the  predictive  process. 

1 .  The  impact  pressure  pulses  have  a  single  wave  form  with  attenuation. 

2.  The  wave  shape  elongates  with  attenuation  so  as  to  conserve  momentum  within 
the  pulse. 

3.  Lateral  relief  processes  determine  the  wave  shape  in  the  thicker  concrete  targets. 

The  range  of  speeds  of  interest  lies  between  the  free-fall  speed  of  the  weapon  and  an 
upper  limit  of  2000  ft/s.  Most  of  the  data  for  the  validation  of  predictions  lie  near  1000  ft/s. 


PREDICTIVE  MODEL 

This  section  describes  a  model  for  characterizing  the  susceptibility  of  a  given 
explosive  to  premature  reaction  during  penetration  through  either  a  concrete  or  a  steel 
barrier.  This  predictive  model  presumes  the  warhead,  limitations,  assumptions,  and 
simplifications  described  above.  Figure  4  is  a  flow  chart  of  the  predictive  model.  After  the 
input  of  data  on  the  barrier,  the  warhead  case,  and  the  explosive,  the  flow  chart  has  two 
parallel  branches  for  the  steel  and  concrete  barriers.  The  need  for  separate  treatments  arises 
because  the  duration  of  the  pressure  pulse  in  the  concrete  thickness  depends  only  on  the 
lateral  relief  processes  while  the  steel  thickness  typically  determines  the  pulse  duration. 
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DATA  INPUT 


FIGURE  4.  Flow  Diagram  for  the  Predictive  Model. 

The  steel  branch  characterizes  the  explosive  by  calculating  a  set  of  critical  impact 
speeds  and  the  corresponding  barrier  thicknesses  for  which  the  given  explosive  can  exhibit 
premature  reaction.  The  prediction  begins  with  the  calculation  of  the  shock  parameters 
generated  by  impact  at  a  given  impact  speed  (starting  with  the  minimum  speed  of  interest). 
The  procedure  then  determines  shock  attenuation  and  transmission  into  the  explosive.  If  the 


NWC  TP  6714 


shock  pressure  lies  below  the  critical  pressure,  Pc,  for  the  SDT  regime  then  the  shock  belongs 
to  the  LALDS  regime  and  the  calculation  of  the  corresponding  critical  pulse  duration  comes 
from  the  empirical  formulation  of  the  underwater  sensitivity  tests.  The  propagation  rate  for 
the  given  impact  speed  and  the  thicknesses  of  barrier  determine  the  barrier  thickness  that 
can  provide  this  pulse  duration.  The  process  continues  with  another  iteration  with  a 
predetermined  incremental  increase  in  impact  speed  until  the  pressure  in  the  explosive 
exceeds  the  value  for  entry  into  the  SDT  regime.  At  this  point,  the  predictive  procedure  has 
completed  the  correlation  of  critical  impact  speed  to  thickness  for  the  given  explosive 
material. 


PREDICTIVE  MODEL 


PENETRATING  WARHEAD  DESIGN 

This  section  presents  a  predictive  model  for  the  reaction  of  an  explosive  to  the 
dynamic  loading  from  a  penetrating  impact  through  either  a  steel  or  a  concrete  barrier.  The 
model  incorporates  the  limitations,  assumptions,  and  simplifications  discussed,  thus  giving 
an  appropriate  form  for  the  characterization  of  the  susceptibility  of  a  given  explosive  to 
premature  reaction  as  a  warhead  component.  The  model  assumes  a  cylindrical  steel  case 
surrounding  a  cylindrical  explosive  charge.  It  further  assumes  strengths  and  dimensions  of 
the  warhead  case  such  that  the  case  avoids  any  significant  deformation  during  the  perforation 
of  the  barrier  (Reference  27).  The  barrier  consists  of  either  a  concrete  slab  or  a  steel  plate 
with  material  behavior  identical  to  that  of  the  warhead  nose.  The  range  of  delivery  speeds 
does  not  exceed  2000  ft/s,  and  the  impacts  occur  at  normal  incidence.  With  these  stipulations 
on  material  behavior,  the  warhead  and  barrier  meet  the  criteria  of  the  preceding  section  for 
the  application  of  a  one-dimensional  anaysis  of  the  early  phases  of  the  transfer  of  internal 
stresses  to  the  explosive  fill. 


EARLY  EVENTS  OF  IMPACT 

Figure  5  shows  the  nose  of  a  cylindrical  warhead  with  a  cylindrical  explosive  charge 
at  the  instant  of  first  contact  with  the  plane  face  of  a  concrete  slab  or  a  steel  plate.  The  region 
of  this  system  near  the  axis  of  symmetry  has  planar  symmetry  so  that  a  simple  x-t  diagram 
can  represent  the  motions  of  its  surfaces  and  the  shock  fronts  that  move  through  it.  Consider 
first  the  target  for  which  the  barrier  is  a  steel  plate.  Figure  6  shows  an  example  with  a  thin 
barrier  typical  of  such  a  target.  The  motions  have  been  exaggerated  for  illustrative  purposes. 
(The  figure  actually  represents  a  far  larger  impact  speed  than  this  model’s  limitations  would 
permit.)  Below  t  =  0,  the  target  plate  approaches  the  warhead  nose.  At  t  =  0,  the  contact 
shown  in  Figure  5  begins.  The  surfaces  in  contact  take  on  a  common  speed.  At  the  same  time, 
two  shocks  move  out  on  either  side  of  the  surface  of  contact.  These  shocks  propagate  the 
changes  of  pressure  and  particle  speed  required  by  the  two  bodies  for  compatible  conditions  at 
the  the  surface  of  contact. 
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The  shock  front  in  the  nose  plate  travels  to  the  interface  with  the  explosive  where  it 
generates  new  reflected  and  transmitted  shocks.  The  transmitted  shock  begins  the  dynamic 
loading  of  the  explosive.  Meanwhile,  the  shock  front  in  the  barrier  travels  to  the  back  where 
it  encounters  a  free  surface  from  which  it  generates  a  reflected  wave  that  reduces  the  pressure 
to  zero  and  establishes  a  particle  speed  twice  that  of  the  incident  wave.  The  reflected  wave 
reaches  the  contact  surface  at  a  time  2htj/Cb,  where  Cb  is  the  propagation  speed  in  the  barrier 
plate  which,  in  general,  depends  on  the  particle  speed.  A  new  transmitted  and  reflected  shock 
should  occur  at  this  interface,  but  if  the  nose  plate  and  the  barrier  have  the  same  material 
properties,  the  shock  will  pass  directy  across  the  interface,  bringing  the  stress-free  state  into 
the  nose  plate.  The  end  of  the  dynamic  loading  of  the  explosive  will  occur  at  an  interval 
2hb/Cp  after  it  began.  Under  these  idealized  plane  shock  loading  conditions,  the  explosive  will 
experience  the  intensity,  P,  of  the  shock  transmitted  into  it  for  a  total  duration  time,  T. 

As  described  already,  this  phase  duration  will  suffer  influences  from  changes  of  wave 
shape  and  internal  processes  of  attenuation  that  accompany  propagation,  but  with 
appropriate  corrections,  the  double  transit  time  provides  a  prediction  of  the  pulse  duration 
that  is  needed  for  the  prediction  of  initiation. 

A  concrete  slab  does  not  have  such  a  simple  early  phase  of  loading  because  it  has  such 
a  great  thickness.  In  a  concrete  barrier,  lateral  relief  effects  arrive  before  the  reflected  wave 
from  the  back  of  the  barrier.  The  time  interval,  D/2Cp,  that  a  relief  wave  of  speed  Cp  requires 
to  reach  the  center  of  the  warhead  of  diameter  D  provides  a  rough  estimate  of  the  pulse 
duration.  Lateral  relief  processes  lack  the  simplicity  of  the  longitudinal  shock  and  its 
reflection. 


CHARACTERIZATION  OF  EXPLOSIVE’S 
SUSCEPTIBILITY  TO  PREMATURE 
REACTION 

The  described  sequence  of  events  applies  to  all  acceptable  variations  of  the  system  and 
impact  speed.  Each  of  these  variations  will  result  in  shock  parameters  P  and  T  that  either 
exceed  or  do  not  exceed  the  critical  conditions  of  Equation  17  and  therefore  result  in  either  a 
premature  reaction  or  no  reaction.  For  any  given  barrier,  warhead,  and  explosive,  the  model 
uses  Equation  17  to  assign  a  single  value  of  impact  speed  that  separates  impacts  with 
premature  reactions  from  those  without  reactions.  The  correlation  of  critical  speeds  to  barrier 
thickness  comprises  a  measure  of  the  susceptibility  of  the  explosive  fill  to  premature  reaction 
that  has  its  formulation  in  terms  of  the  conventional  parameters  for  measuring  the  resistance 
of  the  barrier  to  penetration  and  the  parameter  for  measuring  the  penetration  potential  by 
the  warhead.  This  characterization  applies  to  a  particular  nose  plate  thickness  for  which  it 
provides  the  information  required  to  predict  the  capacity  of  a  particular  explosive  to  survive 
impact  against  steel  plate  barriers.  Such  a  characterization  also  permits  the  comparison  of 
different  explosives  in  the  same  application. 

SIMPLE  ALGORITHM 

The  task  of  predicting  the  minimum  impact  speed  for  premature  reaction  in  a  fixed 
warhead-barrier  system  requires  finding  the  speed  for  which  the  double  transit  time  and  the 
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transmitted  pressure  simultaneously  satisfy  Equation  17.  The  equations  are  readily  written, 
but  it  is  quickly  apparent  that  these  equations  do  not  reduce  to  an  explicit  expression  for 
closure  speed,  V,  in  terms  of  the  other  parameters.  The  solution  requires  some  kind  of 
iterative  process.  The  closely  related  task  of  finding  combinations  of  speed  and  barrier 
thickness  that  satisfy  Equation  17  has  a  far  more  direct  and  simple  procedure.  Given  an 
impact  speed,  the  equations  for  the  generation  of  impact  pressure,  its  propagation  to  the 
metal-explosive  interface,  and  transmission  into  the  explosive  (Equations  6  through  15) 
provide  a  straightforward  means  for  determning  the  pressure  in  the  explosive.  The  equation 
for  the  initiation  of  reaction  (Equation  17)  determines  the  critical  pulse  duration  required  for 
initiation;  the  barrier  plate  thickness  that  provides  this  pulse  duration  comes  from  the 
equation  for  the  double  transit  time,  T  =  2hb/Cfc>. 

The  following  algorithm  used  this  approach  to  determine  a  correlation  of  impact  speed 
to  barrier  thickness  corresponding  to  the  critical  conditions  for  the  initiation  of  reaction  in  the 
explosive. 

1.  Make  a  systematic  choice  of  impact  speed,  V. 


Example:  Choose  1000  ft/s  (300  m/s)  for  the  initial  speed.  If  Px  ^  Pc  and  T  <  1  ms, 
increase  subsequent  speed  by  a  predetermined  amount,  AV  until  T  >  1  ms. 

2.  Determine  the  shock  parameters  at  the  contact  surface. 


u  =  V/2 

(18) 

n 

P  =  p  (a  +  6  -  )v/2 

(19) 

n  r  n\  n  n  2  / 

3.  Estimated  attenuated  values  of  shock  parameters  at  the  metal-explosive  interface. 

(20) 


-ah 

P  =P  e  n 

n  n 


u  =  (  —  p  a  +  'Joia~  +  4p  b  P  )/  2p  b 

n  \  r  n  n  r  n  n  r  n  n  n  J'  n  n 

4.  Calculate  shock  parameters  transmitted  across  the  metal-explosive  interface. 

|-B+  Vff2  -  4AC 


'21) 


(22) 


u  = 

x 


2A 


where 

A 

B 

C 
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and 


P  =  (a  +  b  u  )u  p 

X  X  X  X  X  X 

5.  Determine  the  critical  pulse  duration. 


P" 

X 


(23) 


(24) 


6.  Calculate  the  barrier  thickness  corresponding  to  the  critical  pulse  duration. 


f2ab  +  bbV\ 


(25) 


Repeat  steps  1  through  6  until  both  the  critical  pressure  for  the  SDT  regime  and  the 
pulse  duration  have  been  reached.  The  parameters  and  the  symbols  used  for  them  appear  in 
Table  1. 


Table  1  lists  the  13  parameters  that  characterize  an  impact  system  in  this  predictive 
model.  At  most,  12  of  these  enter  into  any  given  calculation  since  the  barrier  and  the  nose 
plates  have  different  Hugoniot  parameters  only  in  the  case  of  the  concrete  barrier  for  which 
the  thickness,  hb,  does  not  enter  the  calculations.  The  parameters  determine  the  Hugoniots  of 
the  three  materials,  the  thickness  parameters  of  the  one-space-dimensional  approach,  the 
wave  attenuation  exponent,  and  the  explosive  sensitivity  parameters. 


TABLE  1 .  Parameters  of  the  Algorithm  for  the 
Penetration  of  Steel  Plates. 


Parameter  type 

Material 

Barrier 

Nose  plate 

Explosive 

Hugoniot 

a„ 

a. 

bb 

bn 

b. 

Pb 

Pn 

P. 

Attenuation 

Q 

Thickness 

hb 

bn 

Sensitivity 

k 

n 

p 

Shock  Generation 

The  generation  of  plane  shocks  by  the  barrier  and  the  nose  plate  accomplishes  the 
simultaneous  and  instantaneous  accommodation  of  these  bodies  to  the  impact  speed,  V  In 
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the  general  case  of  different  barrier  and  nose  plate  materials,  the  pressure  in  both  waves  has 
a  magnitude  given  by  Equation  5 


P  =  p  (a  +  b  u  )u 

n  r  n  n  n  n  n 


P.  =  p ,(a  +  b  u,)u, 

b  rb  b  n  b  b 


(26) 

(27) 


and  for  compatible  surface  speeds 


u=  V  -  u 
b  n 


(28) 


so  that 


Pb  =  =  h 


ct.+bAV  —  u  ) 

V  -  U 

b  b  n 

n 

(29) 


thus 


(p  6  -  p  b  )u“  -  (2 p  b  V  +  p  a  +  p ,a  )u  +  (p,a,  +  p.b,  V)V  =  0 

rb  b  n  n  n  bo  rn  n  rb  b  n  rb  b  6  6 

and  the  solution,  which  Figure  1  illustrates,  has  the  analytical  form 


(30) 


B  +  VB2-  4 AC  <31) 

u  =  - 

n  2  A 

where 

A  —  pbbb  -  pnbn 
B  =  2pbbbV  +  pnan  +  pbab 
C  =  (pbab  4-  pbbbV)V 

For  the  special  case  of  steel-on-steel  impact,  the  Hugoniots  in  Figure  1  become  mirror  images. 
Equation  30  reduces  by  the  cancelation  of  some  terms  to 


(2p  b  V  +  2p  a  )u  -(p6V  +  pa)V  =  0 


which  has  the  solution  (Equation  18) 


(32) 


u  =  V/2 

n 

so  that  Equation  26  becomes  Equation  19 


P 

n 


=  p  a 
n\  n 


V  \ 

+  b  -  V/2 
n  2  / 
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Changes  of  Wave  Shape  and  Attenuation 

As  the  shock  generated  by  impact  propagates  to  the  metal-explosive  interface,  the 
formation  of  a  multiwave  structure  and  dissipative  processes  both  cause  an  attenuation  of  the 
wave.  (The  present  model  makes  no  attempt  to  account  explictly  for  an  elastic  wave  or  for 
changes  in  the  shape  of  the  high  amplitude  wave  due  to  phase  changes.)  All  attenuation 
effects  appear  as  the  exponential  decay  of  Equation  14.  Thus,  the  pressure  Pn  at  the  barrier- 
nose  plate  interface  decays  according  to 

~  ~ah,i 

P  =  P  e 

n  n 

In  order  to  account  for  the  lengthening  of  the  pulse  that  accompanies  the  attenuation,  the 
pulse  duration  has  a  compensating  elongation 


T 

n 


so  as  to  conserve  the  impulse  of  the  wave 


P  T  =  P  T 

n  n 


(33) 


(34) 


h*. 


r. 

>. 


v, 


Figure  7  illustrates  some  of  the  details  of  this  simplification.  It  shows  the  principle 
features  of  Figure  6  but  with  the  refinements  of  an  elastic  precursor  and  the  representation  of 
the  wave  reflected  from  the  free  surface  of  the  barrier  by  a  set  of  divergent  lines.  Each  line 
represents  some  particular  relief  pressure  value  that  propagates  according  to  its  intensity. 


FIGURE  7.  An  x-t  Diagram  Illustrating  the 
Simplified  Treatment  of  Attenuation  and 
Change  of  Wave  Shape. 
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The  fan  actually  has  a  continuous  smear  of  wavelets.  Such  a  pictorial  device  more  accurately 
represents  the  rarefaction  process  than  a  single  sharp-fronted  shock.  It  implies  that  the  tail  of 
the  pressure  pulse  spreads  out  in  a  continuous  decline  of  some  sort.  The  present  model 
replaces  these  details  of  rarefaction,  the  effect  of  the  precursor  wave  and  any  other  changes  of 
shape  with  a  single  decay  of  a  single  shock  front  and  a  compensating  pulse-elongation  factor. 
This  simplification  accounts  for  attenuation  effects  and  the  known  tendency  for  pulse 
elongation  to  accompany  attenuation  in  a  simple  way  and  avoids  the  procedural  complexity 
that  a  more  detailed  treatment  would  necessarily  require  (and  then  only  by  using 
considerable  speculation).  The  simplification  does  place  a  heavy  burden  on  the  constant  a. 

Equation  5  for  the  attenuated  shock  becomes 


P  =  p  a  u  +■  p  b  77“ 

n  n  n  n  ri  n  ri 


(35) 


The  solution  of  this  quadratic  equation  in  tTn  gives  the  particle  speed  as  a  function  of  the 
pressure,  Pn,  and  the  Ilugoniot  parameters  (Equation  21 ) 


■pa  + 
n  n 


s/77 


pa  ■+■  4p  b  P 

n  n  ’inn 


"  2p  b 

r  n  n 

This  amounts  to  finding  the  new  particle  speed  on  the  P,u  form  of  the  Ilugoniot. 


Wave  Reflection  and  Transmission 


The  shock  transmitted  into  the  explosive  has  pressure,  1\,  and  particle  speed,  ux, 
obtained  by  the  impedance  matching  process  illustrated  in  Figure  8.  The  reflected  wave 
proves  of  no  consequence  since,  as  can  be  seen  in  Figure  6,  it  has  no  effect  until  well  after  the 
cessation  of  pressure  in  the  explosive  by  the  wave  coming  from  the  free  surface  of  the  barrier 
Continuity  of  pressure  and  particle  speed  require  that 


p  (a  +■  b  u  )u  =  p 

x  x  x  x  x  n 


a  +  b  (2(7  —  u  ) 

n  ri  n  x 


(277  -  u  ) 

n  x 


which  reduces  to 


(p  b  —oh  )u~  +  (pa  -fpa  +  4p  6  77  )u  —  (2p  a  77  +  4p  b  77“)  =  0  (36) 

.tjc  r  n  n  x  r  x  i  n  n  r  n  n  n  x  n  n  n  r  n  n  n 

and  has  the  solution  given  by  Equation  22 


(  -  B  +  v/fr  +  4 AC 


2  A 


22 
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PARTICLE  VELOCITY  (u) 


FIGURE  8.  Hugoniot  Curves  for  Determining  the  Particle  Speed  and 
Pressure  Transmitted  into  the  Explosive  Fill. 


for 

A  =  Pxb*  -  Pnbn 
B  =  Pxax  p„an  "h  4pnbnun 
C  =  2pnanun  +  4pnbnun2 

so  that  the  pressure  in  the  explosive  becomes  Equation  23 

P  =  (a  +  b  u  )u  p 

X  X  X  X  Xr  X 


Critical  Pulse  Duration  for 
Given  Pressure 

The  minimum  pulse  duration  for  sustained  chemical  reaction  under  the  given 
pressure  comes  from  Equation  17  that  summarizes  results  from  underwater  sensitivity  tests. 

k 

T  =  — 

pn 

X 

The  model  treats  this  parameter  as  a  sharply  defined  critical  value  as  a  matter  of  expedience 
It  gives  a  deterministic  result  rather  than  complicate  the  procedure  with  poorly  known 
variance  parameters 
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where  the  mean  propagation  rate  has  the  value  of  the  small  amplitude  longitudinal  wave 
speed  in  steel.  The  following  algorithm  starts  with  this  estimate  of  pulse  duration  and  back- 
calculates  the  impact  speed  that  will  generate  the  critical  pressure  corresponding  to  this  pulse 
duration. 

1.  Calculate  the  relief  time  by  means  of  the  elastic  small  amplitude  propagation  rate 


2.  Calculate  the  pressure,  Px,  for  this  pulse  duration,  T,  and  the  corresponding 
particle  speed, u  x 


—  p  a  +  \/p~a2  +  4p  b  P  ) 

rx  x  rx  X  I  x  x) 


2p  b 

X  X 


3.  Calculate  the  pressure,  Pn,  in  the  nose  plate  at  the  steel  explosive  interface  and 
the  corresponding  particle  speed,  un 


B-\/b2-4AC 


A  =  -4pnbn 

B  =  2pnan  -  4pnbn  ux 

C  =  (Pnan  4"  Pxax)ux  +  (Px^x  -  Pn^nlUx^ 


P  =  a u  +6F  p 

n  \  n  rx  n  n  /  n 


4.  Calculate  the  initial  shock  pressure,  Pn,  and  the  corresponding  particle  speed,  un 


~  ahn 

P  =  P  e  n 

n  n 


(  — p  a  +  \/p “a2  +  4p  b  P  ) 

\  r  n  n  r  n  n  n  n  n  / 


2p  b 

r  n  n 
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5.  Calculate  the  impact  speed,  V,  that  will  generate  the  pressure,  Pn,  and  particle 
speed, un 


(  -B  +  \/tf2-4Ac) 

v  =  - - L 

2  A 


A  =  pbbb 
®  —  Pbab  — 2pbbbun 

c  =  -(puan  +  Pbab)un  +  (Pbhb-Pnbn>un2 


SAMPLE  CALCULATIONS  FOR  THREE  EXPLOSIVES 

Sample  calculations  for  three  explosive  fills  of  a  penetrating  warhead  illustrate  the 
use  of  the  preceding  algorithms.  These  provide  examples  of  the  impact  data,  the  final  results, 
and  how  the  data  and  results  are  presented.  The  example  illustrates  the  variation  in 
performance  of  explosives.  These  sample  calculations  also  indicate  the  kind  of  data  that  are 
needed  to  characterize  additional  explosives 

The  general  features  of  the  penetrator  design  were  described  in  the  Penetrating 
Warhead  Design  section.  The  warhead  and  explosive  parameters  required  for  the  algorithm 
for  steel  plates  and  concrete  barriers  are  listed  in  Table  1  Specific  values  of  the  warhead  and 
plate  and  barrier  parameters  that  have  been  used  in  the  sample  calculations  are  listed  in 
Table  2.  These  specific  values  come  from  the  selection  of  mild  steel  for  both  the  penetrator 
case  and  the  target  plates.  The  penetrator  dimensions  are  those  of  the  heavy  wall  penetrator 
(Reference  28).  The  Hugoniot  parameters  for  mild  steel  are  given  in  Reference  29.  The 
concrete  data  are  taken  from  measurements  obtained  at  NWC  and  reported  in  informal 
reports.  The  attenuation  parameter  is  given  as  a  range  of  values,  0.01  0.03,  estimated  from 
observed  shock  attenuations  in  mild  steel  (Reference  20) 

The  values  of  the  explosive  parameters  used  in  the  sample  calculations  are  listed  in 
Table  3.  The  Hugoniot  properties  of  the  explosives  were  determined  by  wedge  test  data 
reported  in  Reference  22  (H-6)  and  Reference  30  (PBXC-1 17(E)  and  PBXW-109(E)).  The 
values  of  the  sensitivity  parameters  came  from  underwater  sensitivity  test  data  (Reference 
23)  that  have  been  fitted  to  the  form  of  Equation  17  and  are  based  on  extremely  limited  data 
(Reference  31). 
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TABLE  2.  Warhead  Case  and  Barrier 
Parameters. 


Parameter  type 

Material 

Steel 

Concrete 

Hugoniot 

p  =  7.85  g/cm3 

p  =  2.5 

a  =  3.85  km/s 

a  =  2.27 

b  =  1.67 

b  =  2.52 

Attenuation 

a  =  0.01-0.03  (1/mm) 

Thickness 

h  =  38.1  mm 

h  =  609.6  mm 

TABLE  3.  Explosive  Parameters. 


Parameter  type 

Explosive 

H-6 

PBXW-109(E) 

PBXC-1 17(E) 

Hugoniot 

P 

1.71  g/cm3 

1 .66  g/cm3 

1 .77  g/cm3 

a 

1.9  km/s 

1.75  km/s 

2.4  km/s 

b 

1.7 

2.78 

2.47 

Sensitivity 

k 

8.0 

7.19 

6.57 

n 

1.75 

2.46 

2.05 

Pc 

13  GPa 

13  GPa 

13  GPa 

The  parameters  of  Tables  2  and  3  were  used  in  the  algorithms  for  steel  plates  and 
concrete  slabs  as  indicated  by  the  flow  diagram  of  Figure  4.  The  steps  of  the  algorithms  were 
formulated  as  a  program  PEMEX  for  a  Hewlett  Packard  9845  desktop  calculator,  which  was 
used  to  expedite  the  calculations.  A  printout  of  PEMEX  is  given  in  the  appendix.  The  result 
of  these  computations  is  a  locus  of  critical  conditions  of  speed  and  barrier  thickness  for  the 
initiation  of  reaction.  These  are  shown  in  Figures  9  through  11.  These  loci  separate  speed 
and  barrier  thickness  combinations  that  produce  no  reaction,  those  below  the  locus,  from 
combinations  that  produce  an  undesirable  reaction,  those  above  the  locus.  Figures  9  and  10 
also  contain  data  from  heavywall  penetrator  tests  conducted  at  N  WC  (Reference  30) 
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THICKNESS,  mm 


Nose  plate: 

Rhon  =  7.857  g/cm  3 
An  =  3.85  km/s 
Bn  =  1  67 
Thn  =  38.1  mm 


Barrier: 

Rhob  =  7.857  g/cm3 
Ab  =  3.85  km/s 
Bb  =  1.67 


Explosive:  H-6 
Rhox  =1.71  g/cm3 
Ax  =  2.83  km/s 
Bx  =  1.7 

Pc  =  1 3  GPa 
Ks  =  8 
Ns  =  1.8 
Tmax  =  1000  ps 


FIGURE  9.  Critical  Speed  and  Barrier  Thickness  for  the  Explosive  H-6. 
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Nose  plate: 

Rhon  =  7.857  g/cm  3 
An  =  3.85  km/s 
Bn  =  1.67 
Thn  =  38.1  mm 


Barrier: 

Rhob  =  7.857  g/cm3 
Ab  =  3  85  km/s 
Bb  =  1.67 


Explosive:  PBXW- 109(E) 
Rhox  =  1.66g/cm3 
Ax  =  1 .75  km/s 
Bx  =  2.78 


Pc  =  13  GPa 
Ks  =  7  2 
Ns  =  2.5 
Tmax  =  1000  ps 


FIGURE  10.  Critical  Speed  and  Barrier  Thickness  for  the  Explosive 
PBXW*109(E). 


0.5 


Nose  plate: 

Rhon  =  7  857g/cm3 
An  =  3.85  km/s 
Bn  =  1  67 
Thn  =  38  1  mm 


Barrier: 

Rhob  =  7  857  g/cm3 
Ab  =  3.85  km/s 
Bb  =  1.67 


Explosive:  PBXC-1 17(E) 
Rhox  =  1  77  g/cm^ 

Ax  =  2  4  km/s 
Bx  =  2  47 

Pc  =  13  GPa 
Ks  =  6  6 
Ns  =  2 

Tmax  =  1000  ys 


FIGURE  1 1 .  Critical  Speed  and  Barrier  Thickness  for  the  Explosive 
PBXC-1 17(E). 
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CONCLUSIONS  AND  RECOMMENDATIONS 

A  model  has  been  developed  for  predicting  premature  reactions  in  the  explosive  fills  of 
a  standard  warhead  design  during  perforating  impacts.  The  model  requires  the  equations  of 
state  of  warhead  structural  materials  and  data  on  the  explosive  fills  derived  from  underwater 
sensitivity  tests  and  wedge  tests.  The  model  has  been  applied  to  steel  as  the  warhead’s 
structural  material  and  three  explosives  for  which  the  required  data  are  available.  A  limited 
amount  of  full-scale  test  data  has  been  obtained  from  sled  tests  at  SNORT.  These  data  show 
that  the  model  is  consistent  with  this  limited  amount  of  data.  It  is  obvious  that  the  model 
depends  strongly  on  the  effective  attenuation  of  shock  waves  in  steel.  The  attenuation, 
described  by  an  exponential  decay,  is  actually  a  complex  phenomenon  based  on  the 
elastic/plastic  behavior  of  steel.  From  Figures  9  and  10  we  can  see  that  a  range  of  values  of 
the  attenuation  factor  are  consistent  with  the  model.  As  more  experiments  are  done,  it  may 
be  possible  to  determine  the  effective  attenuation  factor  more  precisely. 

The  model  provides  a  means  for  comparing  various  explosive  fills.  It  also  provides  a 
correlation  to  laboratory  tests  and,  therefore,  affords  a  means  for  screening  explosives  for  full- 
scale  testing  and  the  data  from  laboratory  testing.  Clearly,  there  is  a  need  for  further  full- 
scale  testing  in  which  the  model  can  assist  in  developing  a  data  base  on  the  susceptibility  of 
explosives  to  premature  reaction  during  penetrating  impacts.  Through  such  a  testing 
program  the  model  can  be  more  fully  verified  and  modified 

The  model  can  also  be  applied  to  theoretical  analysis  using  finite  element  or  finite 
differences  modes  of  numerical  analyses.  Analyses  that  have  already  been  performed  have 
established  the  feasibility  of  this  approach.  One  of  the  features  revealed  by  the  analyses  is 
deformations  of  the  explosive  fill  that  results  in  separation  from  the  rear  surface  of  the 
warhead  followed  by  recovery  with  significant  impact  speeds  on  recovery.  The  model  provides 
a  means  for  assessing  the  potential  for  the  initiation  of  violent  reactions  in  this  process 
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Appendix 

PRKMKX COMIMJTKK  OUTPUT 


PREMEX  convert*  Underwater  Seninivity  Test  Data  arid  Wedge  Test 
into  predictions  of  premature  explosive  reactions  that  develop  1 
penetrating  impacts  against  a  variety  of  barriers  by  flat-ended 
warheads  striking  at  normal  incidence. 


DIMENSIONED  VARIABLES 


DIM  R*C 10] , Vg<201 > , Tbgf 201 > , Ed* I  121 . Nd*C 121 , Bd*C 12] , U*C*J 


OUTPUT  OPTIONS 


Menu:  PRINTER  IS  16 
PRINT  PAGE 

* 

PI ag= 1 
Opt  =0 

PRINT  "OPTION  FUNCTION" 

PRINT  L  I  N  <  1  > 

PRINT  "  1  Correlation  of  Critical  Speeds  to  Critical  Barrier" 

PRINT  "  Thicknesses  for  Thin  Metal  Barriers" 

IF  0pt=l  THEN  Begin 
PRINT  LIN'  1 

PRINT  "  2  A  Single  Critical  Speed  for  a  Thick  Slab  Barrier" 

IF  Opt =2  THEN  Begin 
PRINT  L I N  < 1  • 

PRINT  "3  A  Single  Critical  Speed  for  a  Thick-metal  Barrier" 

PRINT  “  of  Given  Thickness" 

IF  0pt=3  THEN  Begin 
PRINT  LINC3) 

INPUT  "OPT  I  ON  ",  Opt 

Opt  =  I  NT c Opt  >  MOD  4 
IF  Opt -0  THEN  GOTO  Menu 
PRINT  PAGE 

PRINT  "OPTION  FUNCTION" 

ON  Opt  GOTO  210,250,280 

1  BASIC  INPUT  FOR  PROBLEM 


Begin:  1 

No  =  e  _p  1  at  e  :  PRINT  L  I N 1  2  » ;  "  Nose-  p  late  c  h  ar  ac  t  er  i  s i  c  s  " 

INPUT  "Nose-plate  designation  rl2  c h ar ac ’ er s 1 " , Nd* 

PRINT  “Nose-plate:  ";Nd* 

INPUT  "Nose-plate  density  ig«i/cc)",Rhon 
PRINT  "Phon=  "  ;  Rhon;  "  i  ginrc  ) " 

INPUT  "Nose-plate  Hugomo*  A,B  (  Km  sec  ,  -  1 " ,  An,  Bn 
PRINT  " An=  "  ;  An ;  "  »  k'rii  sec  '  B ri  =  "  j  Bn 
INPUT  “Nose-pla'e  thickness  ( mm i " , Thn 
PRINT  "Thn=" ; Thn; ”  i mm > " 

IF  0pt=2  THEN  INPUT  "Nose-plate  diameter  'mm.'",Dia 
IF  Opt  *2  THEN  PRINT  "Dia="i Di a; “  tmi" 

IF  (Opt =2'  AND  '  D  l  a=0  '  THEN  GOTO  520 
INPUT  "Nose-plate  attenuation  'l  mm'". Ain 
PRINT  "A  1  n=  "  ;  A  1  n  ;  "  1  !■  m  n,  I  " 

IF  0pt=2  THEN  INPUT  "Nose-plate  sound  speed  'km  sec '",Cp 
IF  Opt  *  2  THEN  PRINT  "Cp=";Cp;"  (Km  'sect" 

IF  <0pt=2'  AND  <Cp=0t  THEN  GOTO  570 
Barrier:  IF  Opt ; 1  THEN  GOTO  Errr 
Rhob=Rhon 
Ab  =  An 
B  b  *  B  n 
B  d  *  -  N  d  * 

Brrr  :  1 

R*  =  " Bar  r i er " 


>  .  •  .*•  .N  , 


‘,1  •/  V  V  V  's  "ji  *  •  •»'  ,t  jN  ,% 


■  v'vV  .  S- - 

.1  m'>  ck  V 
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670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 


IF  Opt  =2  THEN  R«="S1 ab" 

PRINT  LIN<2);R#; "  characteristic 
IF  Opt>l  THEN  DISP  "Barrier 
THEN  INPUT  Bd* 


desigr.it  ion  <12  characters)"; 


IF  Opt>l 
PRINT 
IF  Opt>l 
PRINT 
IF  Opt>l 
PRINT 


"Barrier:  “;Bd* 

THEN  INPUT  "Barrier  density  c gm  c c > " , Rhob 
“  Rhob= " ; Rhob ; “  cgmrcc)" 

THEN  INPUT  "Barrier  Hugomots  P,B  (Km-sec,- 
"Pb*” ; Pb; " <Km/sec >  Bb=“;Bb 
IF  0pt=3  THEN  INPUT  "Barrier  thickness  (mm>",Thb 
IF  Opt  =3  THEN  PRINT  "  Thb  =  “ ; Thb ; "  <mm>“ 

IF  < Opt  =3 )  PND  < Thb  =  0 )  THEN  GOTO  760 


> " , Pb, Bb 


790  Explosive: 


800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 


PRINT  LIN<2>; "Explosive  characteristics" 

INPUT  "Explosive  designation  <12  c har ac t ers ) " , Ed# 
PRINT  "Explosive:  ";Ed* 

INPUT  "Explosive's  density  < gmx c c > “ , Rhox 
PRINT  "Rhox=“ ; Rhox; "  (g#/cc)“ 

INPUT  "Explosive  Hugoniot  P,B  < Km/sec , - V , Rx , Bx 
PRINT  "Px  =  ";Px;  "(Kin/sec  )  Bx=";Bx; 

PRINT  TPB<40> ;  “  CP-'Ns*T  =  ks,  <P  in  GPP,  T  in  usee)]" 
INPUT  "Critical  pressure  <GPP)",Pc 
PRINT  "Pc  =  ”;Pc; “  ( GPP  > " 

INPUT  "Sensitivity  coeff icient",Ks 
PRINT  "Ks=";Ks 

INPUT  "Sensitivity  exponent", Ns 
PRINT  “ N  s  =  “ ; N  s 

IF  <Opt>l)  PND  <Ns=0)  THEN  920 

BEGIN  THE  CPLCULPT I ONS 
INPUT  "Make  corrections,  then  <CONT>”,R# 


985  Cax=Rhox*Rx 
990  Cbx=Rhox*Bx 

1000  Can“Rhon*Pn 

1010  Cbn=Rhon*Bn 

1020  Cab=Rhob*Pb 

1030  Cbb=Rhob*Bb 

1040 
1050 
1060 

1070  ON  Opt  GOTO  Opt  1 , 0pt2, 0pt3 
1080 
1090 
1100 

1110  Opt  1 
1120 


OPTION  BRPNCH  POINT 


OPTION  1 


1  130 
1  140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 

1230  Opt 1_1 oop: V=V-De 1 v 


Maximum  time  of  interest  (usee ) ” , Tma> 
IF  Tmax<=0  THEN  Tmax=1000 


! 

INPUT 
Px  =  Pc 

GOSUB  Inurse_trnsfr 

PRINT  1_IN<1>;  "  V  =  "  ;  V ;  "  Kmrsec" 

INPUT  "P  better  value  for  V’",V 

Del v=Vx200 

V=V+Del v 

MPT  Vg=<0) 

MPT  Tbg=  <  0 ) 

Mk»l 


1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 


IF  v<=0  THEN  GOTO  Drau_graph 
Un=V/2 

GOSUB  D i rec t_t rnsf r 

IF  T i me >Tmaxx 1 E6  THEN  GOTO  Draw_graph 
GOSUB  Doub 1 e_t r vs  a  I 
Vg<  Mk ) = V 
Tbg<Mk )*Tb 

IF  Mk=201  THEN  GOTO  Dr au_gr aph 
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.N  .v  . 
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Mk  =Mk ♦  1 

GOTO  Opt  1  1  oop 


OPTION  2 


T  i  me  =  t)  l  a/(2*Cp)*lE-6 

Px  =  <Ks-'<Ti  me*lE6>X<:  1/Ns> 

GOSUB  I  nvrse_t rnsfr 

U*  = " Unref i ned" 

GOTO  Print_data 
0pt2_back:  Ut=" Ref i ned" 

Loops=0 

0pt2_loop:  GOSUB  Bar>_impmat 

GOSUB  Di rec t_t rnsf r 
GOSUB  Variab_rel lef 

IF  RBSCDia-DxXBiaxlES  THEN  GOTO  Prir,t_data 
V«V+V*(Dx-Di  aX(Ns*Di  4)x2 
Loops=Loops+l 

IF  Loops<200  THEN  GOTO  0pt2_!oop 
PRINT  ”V=";V 

IF  Loops<250  THEN  GOTO  0pt2_loop 
PRINTER  IS  0 

PRINT  L I  H<  1 ) ;  “Ref  i  n«d  critical  speed  did  riot  converge 


PRINTER  IS  16 
GOTO  Menu 

OPTION  3 


INPUT  "What  is  your  estimate  for  V,,‘,V 
Loops=0 

0pt3_loop:  GOSUB  Ban_impmat 

GOSUB  D i rect_trnsf r 
GOSUB  Doub 1 e_t rusal 

IF  RBS<Thb-TbXThbxlE9  THEN  GOTO  Print_data 
V»  V  + V*  <  Tb-Thb  >  -'  (  Ns*Thb  >  x2 
Loops=Loops+ 1 

IF  Loops<200  THEN  GOTO  0pt3_loop 
PRINT  "V=" ; V 

IF  Loops<250  THEN  GOTO  0pt3_loop 
PRINTER  IS  0 

PRINT  L I N ( 1 ) ; "Cr i t i c a  1  speed  did  not  converge. . .  ";LI 


PRINTER  IS  16 
GOTO  Menu 


INVERSE  SHOCK  TRRNS VERSE 


I nur se_t r nsf r :  1 

R  =  Cbx 
B  =  Cax 
C  =  -Px 

IF  fl  =  0  THEN  U  x  =  -  C  B 
IF  fl  =  0  THEN  1870 


1860  IJx  =  <  -  B+SQR<  B  “  2-4»R*C  >  X  R,  2 

1870  fl=4*Cbn 

1880  B=2*Can-4»Cbn*Ux 

1390  C  =  -<.  <  Can  +  C  ax  '*Ux+(-Cbn  +  Cbx  )*Ux'2> 

1900  IF  R*0  THEN  Unt  =  -C  B 

1910  IF  R«0  THEN  1930 

1920  Unt  =  < -B+SQR  ’  B~2-4*fl*C  •  R  '2 

1930  Pnt  =Cari*Unt  +Cbn*Unt  2 

1940  Pn«Pnt  *EXP<  R1  n*Thn 

1930  R«Cbn 


N'WC  TP  67 14 


I960 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 
24  30 
2440 
2450 
2460 
2470 
2480 
2490 
2500 
2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 


B  =  Can 
C  =  -Pn 

IF  8=0  THEN  Un  =  -OB 
IF  8=0  THEN  2010 
Un=(  -B  +  SQR<8-'2-4*0*C>  >x8-'2 
IF  Opt  =  1  THEN  V  =  2*Un 
IF  Opt  =  1  THEN  RETURN 
8=Cbb 

B=Cab-2*Cbb*Un 

C  =  -<  <Cab  +  Can)*Un  +  <Cbri-Cbb)*Un^2> 

IF  8=0  THEN  V=-C/B 
IF  8=0  THEN  2090 
V  =  <-B  +  SQR<B''2-4*8*C>  >/8/2 
RETURN 

►  ♦♦♦•in 

B8RR I ER-NOSE  IMPED8NCE  M8TCH 


Ban_i mpmat :  * 

8=Cbb-Cbn 

8=-<Cab+Can+2*Cbb*V) 
C=Cab*V+Cbb*V~2 
IF  8=0  THEN  Un=-CxB 
IF  8=0  THEN  2200 
Un=<-B-SQR<B"2-4*8*C)). I 
RETURN 


/  8 


DIRECT-SHOCK  TR8NSFER 


D i rec  t_trnsfr!  1 

Pn  =  Can*Un*Cbn»Uri-'2 
Pnt=Pn*EXP<-81n*Thn> 

8  =  Cbn 
B  =  Can 
C«-Pm 

IF  8=0  THEN  Unt=-C'B 

IF  8=0  THEN  2340 

Um  =  <  -  B*SQR<B'2-4*8»C))  '8-'2 
8=Cbn-Cbx 

B  =  -  c  Carv»Cax +  4  *Cbn*Unt  ) 

C  =  2*Can  +  Urit*4*Cbn*Unt''i 
IF  8=0  THEN  Ux  =  -C  B 

IF  8=0  THEN  2400 

U>  =  <  -  B-SGR'  B  "  2-4*8»C  >  '8  2 
Px=Cax*Ux+Cb'»Ux~2 
T i »«=Ks  Px " Ns  "  1 E6 
RETURN 


i 


►  +■*♦*  -4 


DOUBLE  TRAVERSAL 

Doub 1 e_t  rv sal :  1 

T b  =  T  i  m e ♦  <  <  A n ♦  B n * U n  '  X’  * E X P 1  - A  1  n  +  Thn  1  ♦  l  E 6 
RETURN 

►♦♦♦♦♦♦♦*4+****»*4*»4»**+*X4*t*t****+****^ 

VRRI8BLE  RELIEF  R8TE 

Var i ab_re  1  i  eF :  1 

Dx  =  2*  T  1  n.e  ♦  <>  8n+  <  Bn-  1  ♦  Un*  »'8n»2*Brv*Un^  '  Rn  +  Bn*Un  '  *  1 E6 
RETURN 

!  CR8PHICS  OUTPUT 


Draw_graph:  1 

IF  Keep=t  THEN  GOTO  L'mits_;e' 
PLOTTER  IS  13, "GRRPHICS" 

Vmax  =  . 0000000000  I 
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2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2755 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2330 
2840 
2850 
2860 
2870 
2880 
2890 
2900 
2910 
2920 
2930 
2940 
2950 
2960 
2970 
2990 
2990 
3000 
3010 
3020 
30  30 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 

r :  “ ; 
3150 
r  l  er 
3160 
"  ‘  u 
3170 
3190 
3  1  90 

3200 
3  2  1  0 
m  c  r 


IF  MflXt  VgtMk -1  )  ,  Vg>  1  )  ><=Vmax  THEN  2650 
Vrnax  = 1 O* Vrnax 
GOTO  2620 

IF  Vrnax  '2- MflXf  Vg>  Mk  -  1  '>  ,  Vgi  1  >  >  THEN  2630 
Vmax  = Vmax/2 

GOTO  2650 
Tbm ax =.0000000001 

IF  MfiX<Tbg<Mk-l>, Tbg< 1 > ><=Tbmax  THEN  2720 
Tbmax=10*Tbmax 
GOTO  2690 

IF  Tbmax-'’  2<  M  AX  ( Tbg  <  Mk  - 1  t ,  Tbg  <  1  >  >  THEN  2750 
T  b  m  ax  =  T  b  m  ax  •-  2 
GOTO  2720 
Limit  s_set :  1 

Tbma*=250 

GRAPHICS 

LOCATE  0,120,0,100 

SCALE  - . 2* Tbmax , 1 . 1 tTbmax , - . 2* Vrnax , 1 . 1  * Vm ax 

LINE  TYPE  3 

CLIP  0 .  Tbmax ,  0 ,  Vrnax 

IF  K  e  e  p  <  >  1  THEN  GRID  Tbmax/ 1 0 , Vrnax/ 1 0 , 0 , 0 
LINE  TYPE  1 

CLIP  2*Tbm  ax  ,1*1  *  Tbmax  f  - .  2*  Vrnax ,  1 .  1  *  Vma.x 

AXES  Tbmax/ 10, Vmax/10, 0,0 
HOVE  Tbg< 1 > , Vg< 1 > 

FOR  L  =  2  TO  Mk - 1  STEP  1 
DRAW  Tbg>L>,VgtL') 

NEXT  L 

MOVE  -Tbmax .'20, -  Vrnax- -20 
CSI2E  3.5 
LABEL  "0" 

MOVE  . 9  7  5  *  T  b  m  ax , - V  m  ax  x 1 0 
LABEL  VAL*  <  Tbmax ) 

MOVE  -T bmax.'  1 0 ,.  99*  Vrnax 
LABEL  VAL*  <  Vrnax  > 

MOVE  . 9* Tbmax '2 , -Vrnax/ 1 0 
LABEL  " Th  ( mm i " 

MOVE  -Tbmax x 1 0 , . 9* Vrnax x 2 
LDIP  PI  2 

LABEL  "V  (Km-'ssc  )  " 

LDIR  0 

MOVE  Tbg< 1 > , 1 . 01+MI1K Vrnax , Vgt 1 ) ) 

IF  K«ep=l  THEN  MOVE  1 . 01 -MIN'  Tbma- ,  Tbg' Mk  -  1  <  ,  Vgv Mk-1  j-' 1 . 02 
CSIZE  2.63 
LABEL  Ed* 

DUMP  GRAPHICS 
EXIT  GRAPHICS 
GOTO  Print  _d at  a 

I  4r-*****JiJmt****-*-*-****^-*-***^*************-************-*****^********** 

1  PRINTED  OUTPUT 

Print_datal  1 

PRINTER  IS  0 

IF  Opt  =  3  THEN  PRINT  L I N <  1  '{"Critical  speed  for  thick-metal  barrie 
1 ;  "  'hfi/sec  >  " 

IF  Opt  =2  THEN  PRINT  LIN'I'JU*;"  critical  spieed  for  thick-slab  bar 
" ; v ; “  ffm-itc  '  “ 

IF  Opt  l  THEN  PRINT  " P  *  « “ ; P  x ; "  ■ GPA  • " 1 T AB <  34  ) ; " T i ms =  ";Tim«*lE6; 

ft  )" 

PRINT  LIN1  1  ' 

IF  'Opt  =2  1  AND  ■ U*  =  "Ref i ned"  '  THEN  GOTO  3  330 

PRINT  ‘'Nose-plats:  "  1  N  d  t ;  T  A  E  ■  25  ■;  "Barrier:  "  ;  Ed* ;  T  AB  ■  50  • ;  "  E  p  1  o  i  i 

;  Ed* 

PRINT  LIN'  1  1  ; 

PRINT  "  Rhori  =  "  ;  R'Hon  1  T  AE  ■  25 ' 1  "  Rhob=  “  ;  R’hob  1  T  AE  '  50  '  1  "  R' Ho  ■  =  "  ;  Rho  "  1  q 
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PRINT  "An=“ ;  fin;  TAB < 25) ;  "Ab="  ;  fib; TAB (50) ;  "  Ax=" ;  fix;  “  < Km-  sec  > 
PRINT  "Bn="  ; Bn; TAB (25); "Bb=" ; Bb; TAB (50); "Bx  =  " ; Bx 
PRINT  "  Thn-  “  J  Thn ; 

IF  Opt  =3  THEN  PRINT  TAB ( 25 > ; " Thb= “ ; Thb ; 

PRINT  “  (in)" 

IF  Opt  =2  THEN  PRINT  “ D i a= “ ; D i a; "  (mm)"; 

PRINT  TflB<50); "Pc  =  “;Pc;  “  ( G P A ) " 

PRINT  "fll  n  =  "  ;  fll  n;  "  <  1  /n )  “  ;  TAB < 50 )  ;  " Ks* "  ;  Ks 
IF  Opt  =2  THEN  PRINT  "Cp=“;Cp;"  (Km-'sec)"; 

PRINT  TfiB( 50 ) ; "Ns= " ; Ns 

IF  Opt  =  1  THEN  PRINT  TAB< 50 ) ; " Tmax  = " ; Tmax ; "  (usec)“ 

PRINT  LIN(l) 

IF  <  Opt  <  >2 )  OR  (U$="Ref  i  ned"  )  THEN  PRINT  LINO) 

PRINTER  IS  16 

IF  < Opt  =2 )  AND  (U*="Unref ined")  THEN  0pt2_back 
PRINT  PAGE 

********************************************************************* 
CONTINUE  CALCULATIONS 

****************************************************************** 
IF  Opt  =  1  THEN  Keep=0 

IF  Opt=l  THEN  INPUT  “Enter  1  to  overlay  next  curve", Keep 
GOTO  Menu 
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NOMENCLATURE 

The  symbols  used  in  the  text  and  equations  of  this  report  are  summarized  below  in 
alphabetical  order. 

Ac,  Bc,  Cc  Terms  in  the  quadratic  formula  for  the  solution  of  Equations  4  and  9 

An,  Bn,  Cn  Terms  in  the  quadratic  formula  for  the  solution  of  Equation  31 

Ax,  Bx,  Cx  Terms  in  the  quadratic  formula  for  the  solution  of  Equation  22 

A,  B,  C  Terms  in  the  quadratic  formula  for  the  solution  of  Equation  45. 

a,  b  Hugoniot  constants 

ai,  bj  Constants  for  the  reflection  Hugoniot 

a2,  bo  Constants  for  the  Hugoniot  in  transmission 

at,,  b|,  Constants  in  the  barrier  material 

an,  bn  Constants  in  the  warhead  nose 

ax>  bx  Constants  in  the  explosive  fill 

bn  Density  in  a  shock  in  the  warhead  nose 

bx  Density  in  a  shock  in  the  explosive  fill 

C  Propagation  speed  through  a  material 

Cb  Propagation  speed  through  a  barrier 

Cp  Longitudinal  sound  speed  in  steel 

D  Diameter  of  warhead 

E  Internal  energy  behind  a  shock  front 

E0  Internal  energy  ahead  of  a  shock  front 

hb  Thickness  of  the  barrier 

hn  Thickness  of  the  warhead  nose 

k  Constant  in  the  relation  between  pressure  and  pulse  duration 

given  by  Equation  17 

n  Constant  in  the  relation  between  pressure  and  pulse  duration 

given  by  Equation  17 

P  Pressure  behind  a  shock  front 

P0  Pressure  ahead  of  a  shock  front 

Px  Pressure  in  a  shock  in  the  explosive 

P  Pressure  in  a  shock  after  attenuation 

Pn  Pressure  in  the  nose  material  of  a  warhead 

Pn  Pressure  in  the  nose  material  of  a  warhead  after  decay  of  the  shock 

Px  Pressure  at  a  shock  wave  in  the  explosive 

T  Pulse  duration 

T2t  Pulse  duration  determined  by  a  relief  wave  that  performs  double  transit  of  a 

given  thickness  of  material 
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Time 

Shock  speed 

Speed  of  the  shock  in  the  material  through  which  an  initial  shock  has  been 
transmitted  and  reflected  from  a  second 
Particle  speed  in  body  1  during  an  impact 

Shock  speed  in  a  material  shocked  by  transmission  from  another  material 
Particle  speed 

Particle  speed  in  body  2  during  an  impact 
Particle  speed  in  the  barrier  during  an  impact 
Particle  speed  in  the  warhead  nose  during  an  impact 
Attenuated  particle  speed  in  the  warhead  nose 
Particle  speed  in  the  explosive  fill 
Impact  speed 

Initial  speed  for  premature  reaction 
Increment  of  impact  speed 
Attenuation  constant 
Density  behind  a  shock  wave 
Density  ahead  of  a  shock  wave 

Density  in  the  shock  developed  in  body  1  during  impact 
Density  in  the  shock  developed  in  body  2  during  impact 
Density  in  a  shock  in  the  barrier 
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